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Abstract 

We consider interactions of fermions with the domain wall bubbles produced during 
a first order phase transition. A new exact solution of the Dirac equations is obtained 
for a wall profile incorporating a position dependent phase factor. The reflection 
coefficients are obtained. 
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1 Introduction 



During the last few years, a considerable amount of work has been dedicated to the possi- 
bility of generating a sizeable baryon asymmetry on the electro-weak phase transition (see 
for example [|, |||, |). For an excess of baryons to develop in an Universe which initially 
has zero baryon number, the already well known following conditions, first enucianted by 
Sakharov, must be met: 1) Some interaction of elementary particles must violate baryon 
number. 2) C and CP must be violated in order that there is not a perfect equality be- 
tween rates of AB ^ processes, since otherwise no asymmetry could evolve from an 
initially symmetric state. 3) A departure from thermal equilibrium must play an essential 
role, since otherwise CPT would assure compensation between processes increasing and 
decreasing the baryon number. Remarkably, the standard model of weak interactions may 
provide all the necessary ingredients for baryogenesis. In particular the third condition 
can be met if the weak phase transition is at least weakly first order. 

In a first order phase transition, the conversion from one phase to the other occurs 
trough nucleation. This happens when the system is either supercooled or superheated. 
Bubbles of the "true" phase (with an expectation value of some Higgs field v ^ 0) expand 
rapidly absorbing the region of the "false" phase (v ^ 0). At the bubble surface there 
is a region or "wall", in principle of microscopic dimensions, which separates the phases. 
The speed of the expanding bubble walls could be in the range 0.1 — 0.9 c Particles in 
the "false" (higher temperature) phase are reflected off the advancing bubble walls, while 
most particles in the low temperature phase are unable to catch up with the receding 
walls, and cannot equilibrate across the phase boundary. Thus one has a departure from 
equilibrium and a baryon asymmetry can be generated. 

In the physical conditions of the early Universe the fermions moving through the 
bubble wall will interact also with the particles in the surrounding plasma, thus a full 
transport problem must be considered. A useful simplifying assumption is to decompose 
the process into two steps, one describing the production of the CP asymmetry on the 
transmission/reflection coefficients when the quarks/antiquarks are scattered on the wall, 
the second describing the transport and the eventual transformation of the CP asymmetry 
into a baryon asymmetry via the baryon number anomaly. 

Assuming that the scattering from the wall is little affected by diffusion corrections, 
the effects of the surrounding plasma can be partially incorporated by introducing a Higgs 
field effective potential taking into account finite temperature corrections to the tree-level 
potential. The structure of the wall depends on this effective potential on a complicated 
way. Fermions passing through the domain wall acquire a mass which is proportional to 
the vacuum expectation value (VEV) of the Higgs field, which is determined from the 
equations of motion of the finite temperature effective action of the bubble. The problem 
of computation of transmission coefficients reduces to the solution of a Dirac equation 
with a space dependent mass term. Exact solutions has been obtained only for two simple 
cases: for a wall profile approximated by a step function ([@, P]) and for an average Higgs 
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field profile of the type (Q): 

= 



l + tanh(i)) (1) 



where the width of the wall is given by 

2^2 ^2 
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and Mh = \/\v/2 is the Higgs mass. This profile is the analytic solution of the Higgs 
field equation for a bubble with its normal along the z-axis and position at z = under 
some simplifying assumptions. There is not variation of the complex phase of the Higgs 
field through the wall. Numerical solutions have been obtained for some more complicated 
profiles which incorporate locally dependent complex phases and which are considered to 
be "reasonable" enough although not neccesarily solution of any equation of motion ( for 
example in [^). 

For a given wall profile, parallel to the x-y plane and normal to the z-axis in its 
own rest frame, in order to compute the reflection coefficient one need only the plane 
wave solutions of the Dirac equation for particles moving along the z-axis. For any other 
incoming direction, the problem can be reduced to the latter performing an appropriate 
Lorentz boost. Following ||^, 0, it is advantageous to work in the chiral basis, reordering 
the spinorial components the Dirac operator can be factorized into 2x2 blocks. For 
solutions with positive energy E 



iEt I i'liz) \ _ „-iEt 
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(2) 



where ipi and ^^2 are eigenspinors of the chirality operator , 75, for the eigenvalue +1 and 
■03, ^^4 we obtain the two equations 

{id, + Q{z))^i = (3) 
{id, + Q*{z))iJii = (4) 

With 

We will deal in this work with the particular case given by the function 

( \ —\ "iQ exp i(— A^^Az + ^0) exp — Az if 2; > , ^ 

I rriQ exp iOq if 2; < 

All the constants are supposed real. A > 0. This function represents a linear phase 
variation with a global difference of over a distance of the order of the wall thickness 
<5 = 1/A. 
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2 Solving the Dirac equation. 



It would be possible to solve Eqs.(^^) with the function given by Eq.(P), reducing them 
to a Bessel-like differential equation. It is possible and advantageous however to use 
perturbation theory to compute the evolution operator of the system. The summation to 
all orders of the perturbation expansion is possible thanks to the special form of Q{z). 
There are two main advantages in doing so: the first one is that the procedure is easily 
generalizable to any number of dimensions (for example to incorporate mixing between 
generations), the second one is that the evolution operator is computed directly and the 
reflection coefficients are easily given in terms of its components. A similar technique 
has been used already in (||8|, |9|) to compute the neutrino oscillation probabilities in solar 
matter. 

The evolution operator of the differential system (^) is given by the path-ordered 
integral 

U{z,zo) = Pexp-i f dzQ{z) (7) 
In this work we are concerned with an operator Q of the form 

g = Q° + v{z) (8) 

With 

/c, a are in general complex. The real part of a is greater than zero. It can be set 
3?(T = 1 without loss of generality. 

Formally, it is possible to solve Eq.(0) by successive iterations: 



C/(z,zo) = C/(°)(z,zo) + E^^"H^,^o); [/(°)(z,^o) = exp(-iQO(z-zo)) (10) 

n=l 

is the well-known integral 

= {-iT dzn... dziU^iz, Zn)V{zn) . . . U\z2, zi)F(zi);70(zi, zo) (11) 
The domain of integration is defined by 

T = Z > Zn > ■ ■ ■ > Zi > Zq. 

Following the same arguments as in one can show that it is enough to compute U 
in the z — > oo limit. The evolution for finite time can be deduced from the expression for 
this limit. Trought elementary manipulations of Eq.(|lT|), we get the elements of C/*^"^ in a 
basis of eigenvectors of Q^: 

< b I I a >= (-i)" exp(-i(Q|;z - Q°zo)) x 

X J2 d"''"exp(zZ„W6fcl + . . . + i2;itffc(n-l)a)^bfel(2;n) • • • Vfc(„_i)„(zi) (12) 

fel,...,fc(n-l) ^ 
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With Wkik2 = Qki ~ Qk2- Qk °f eigenvalues of Q^. 

Due to the dimensionahty of the problem and the special form for V, the summatory 
m Eq.(ll2D either is zero or reduces to only one term depending on whether n is odd or 



even and on the states a, b. For diagonal terms (a = b) the product of V . . .V is always 
zero when n is odd. When n is even there is a single surviving term. For non-diagonal 
terms, the situation is reversed: the only one surviving term appears for n odd. This 
single remaining term is always of the alternating form . . . ^42^1^21 . • •. So, 

< b I [/(") I a >= {-^reM-^{Qtz - Q^o))/!? x | | ^ T I (13) 

[ Vba (- I fc 1)^ " aj^b 

(n) 

All the functions appearing in the integral 1^^ are of exponential type, the following 
equality (|^]) can be applied: 

In{wi, Wn) = / . . . / dXn . ■ ■ dxi CXp ^ 

J zn J zn 



-l)''exp(zo En Wn) 



Wi{wi + W2) . ■ ■ {Wi +W2.--+ Wn) 

(valid if K u;„ < 0, V n ) (14) 



In our case 



wi + W2 + ■.. + Wj = iwbki - L + iwkik2 - L* + + iwf,(j-i)kj - L 
= iwbkj +niL + n2L* 



(15) 



Wbk(j) can take only the values {0, ±25'}. L is a or o"*. One have for the integers 
ni , n2 : ni = n2 or ni = 722 =b 1 . 

For the diagonal terms, a = b, n even (writing s = 2Ei — a): 



\even 

"■aa 



exp zo(n/2)(s + s*) 



s*{s* + s){s* + s + s*)... ((n/2)(s + s*)) 
exp zo(n/2)(s + s*) 

YlU^evenijmis + S*) UU,odd[^* + " + 

_ exp zo{n/2){s + s*) 

~ {s + s*)"/2 (f )!(s + s*)"/2 nri?[sV(5 + s*) + (/ - 1)] 
_ ex.p zo{n/2){s + s*) 

(s + s*)"(f)'[sV(s + s*)](„,/2) 

For non-diagonal terms, taking a = 1, 6 = 2 to simplify the notation: 



(16) 



{n),odd _ (-1) exp zo((ra/2)(s + s*) + s*) 



ab 



s*{s* + s){s* + s + s*) ... ((n/2)(s + s*) + s* 
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-l)exp zo(n/2)(s + s*) 



np2,e.en(j72)(s + S*) U]=1,om{^* + (j " + S*)/2) 

(-l)expzo(ra/2)(s + s*) 

(s + s*)("-i)/2 (^^) !(s + s*)("+i)/2 n[=t^^^^[^V(s + s*) + 1)] 
(-l)expzo(n/2)(g + s*) 

We have used the Pochammer symbol defined by Eq. (p|). 

So, inserting Eqs.(16,17) in Eq.(13) and taking s + s* = —2, zq = 0- 



(17) 



<l\U\l>-- 



0<J 

1+ E Hr(-i^i') 



n/2 



n=2,ei)en 
oo / I 7„ |2 ' 



1+E 



m=l 



e o-f^i I — 



2 ' 4 



(-2)"(f)![-.V2](n/2) 



m![-sV2](„)^ 



(18) 



< 1 I ?7 I 2 >= 



e-^^(-fe) ^ (_i)-(_|fc|2) 



2iA: ^ l\k 



E 



{n~l)/2 



-1) 



(-2)"(V)'[-^V2]((„+i)/2) 



4 y (m-l)![-sV2]M 



e-^^i^oFi(l--;^ 
s* \ 2 4 



and similarly for the other two matrix elements. The matrix U can be written as 
U{z oo, zq) = exp -iHo{z - Zo)Uredizo) 



(19) 



with 0" = A(l + iA9) , k = —mo exp iOq 

imo exp zOq 
1 + {2E- Ae)i 



G 



F G 

G* F* 



^ (l ,^ I mo 

oFi - + (S - 



(20) 



(21) 
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In the case zq ^ 0, F,G would include a factor exp —2zq in its argument. In these formulas 
E, mo are given in units of A, the inverse of the wall thickness. 

To obtain U, evolution operator for the Eq.(Q), one must make the changes k — > 
k* , A9 — > —AO. The matrix become 

U{z^QO,zo) = eyip-iHo{z - zo)Ured{zo) 

Uredizo) = ^*) (22) 



with 



^ imoexp-iA^ /3 |/c|2\ 

F = oi^if 2 + ^^ + ^)'' 4 ) 

See Appendix A for some new formulas for the absolute values of generalized hyper- 
geometric functions which can be obtained from the general properties of U, U. 

Following the same reasoning used in ([^), using the general properties of the evolution 
operator, the evolution for any finite z is given by the matrix 

U{z,zo) = Us{z)-'Us{zo) = U-],{z)e"''^'^'~''>'>Ured{zo) (24) 
In this work, we will make use only of the infinite time limit. 



3 The reflection Coefficient. Results. 

We will follow the same procedure as in (|j|]) for defining the reflection coefficient. In the 
region z ^ oo, the eigenspinors ip2,ip3 correspond to right-moving particles with chirality 
75 = 1 respectively. The solutions V'i)'04 &re identified with left-moving particles 

with chirality — 1. These states are also eigenstates of the hamiltonian in this region 
The momentum eigenstates for z < are obtained diagonalizing the constant matrix 
Q(0) by 

/cosh^, sinhe\ ,i,h2^,^^ = ^=i=_ (25) 
Vsmh^p cosh^py' ^ p ^[E/ \m\Y -I 



with eigenvalues p = ^\/E'^— \ m 

For left-moving particles incident from the symmetric phase, two components coexist 
at z = oo, the incident particle itself and the reflected one by the domain wall. We define 
reflection coefficients R, R as 

ipsioo) = Ripi{oo); ^^2(00) = R^l^A{oo) (26) 
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Imposing the boundary condition that at z < only a left-moving particle with momentum 
p propagates, one gets the expression: 



^(oo) = U{z co,0)u-\p)4^p{0) (27) 
where ipp is a momentum eigenspinor. Then, the reflection coefficient is given by 



R=^^l^. R = (21 (28) 
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or more explicitly: 



U2I - tU22 ^ 2E.^ G* - tF* ( ) 

Uu -tUi2 F-tG ' ^ ^ ' 



and similarly for R. 

We are interested in the asymmetry between the reflection probabilities. The quantity 
^4=1 i?p — I i?pis displayed in figs.(||-^) as a function of the diverse parameters involved 
A0, 9q and for different values of the product particle mass x wall width | m | (5. 

In Fig.(l^) the asymmetry A is plotted as a function of the dimensionless quantity 
E/ \ m \ and some fixed values of A0. For = vr (plot B) the results obtained here 
coincide with the numerical result obtained in (Fig. 2). As the phase difference 
increases, the height of the peaks keeps unaltered but their position moves according to 
the rule 

{E/m)peak =1 I /{2m5). 
The first parameter of the hyper geometric functions in Eqs.(p^23) is real for these peak 



values. As the value of A0 gets higher some resonant effect becomes apparent, the peaks 
get sharper and the asymmetry is only non-negligible around the peaks. This is particu- 
larly evident in the last plot (C). This effect could lead to an essentially monochromatic 
asymmetry in scenarios with highly oscillatory or random phase differences. 

In Fig.(^) we check the influence of the value of the phase ^o- The plot (A) corresponds 
to plot (C) in Fig. (HI). Here we have set = Some new dips appear at lower E/m, 
otherwise the position and width of the peaks remain unaltered. In the plot (B) we present 
the variation of A (peak values at 9q) in a full cycle of 9q. The phase appears only 
in the multiplicative exponential coefficient of G in Formula (^). From this, the soft, 
quasi-sinusoidal behavior observed in the graphic is understood. 

Finally in Fig.(^) we plot A as a function of the parameter A0, for 9q fixed and for 
different values of the product m5 as before. As expected the peaks in the asymmetry 
function appear for values such that 

As conclusion, we were able to solve the Dirac equation with a space dependent complex 
mass term which, although not a solution of the equations of motion, reproduces the 
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expected variation in module and in complex phase of the average Higgs field across 
the wall. From the solution of the Dirac equation the particle-antiparticle transmission 
asymmetry is computed. The analytical results presented here confirm previous numerical 
computations and predict an unexpected behavior for highly oscillatory phase fields. 

We note that with little modifications our method can be used to solve the Dirac 
equation for a purely local dependent phase mass term (A 0,A9X 7^ in Eq.(^) or 
equivalently a purely imaginary in Eq.(P)). In order to circumvent convergence problems 
this must be done taking the limit — > in the finite propagation time Eq. p^) . 



A Appendix: some old and new formulas about Hypergeo- 
metric Functions 



The generalized Hypergeometric function (||10|) is defined by 

TP ^ , r . ^ (Ql)(n)---(«p)(n) Z'^ 

Ma„ ...,a,M,..., b,; z) = ^ (,^)^^^ . . . (30) 
where the Pochammer symbol is 

=r(n + z)/r(z) (31) 

in particular 

The derivative of qFi is again a 0-^1 function: 

= -oi^il+T;^; (33) 

dz 7 

The function qFi is related to the Bessel Functions by the formula 

For a Q as given by Eq.(^) is traceless, det Q{z,zq) = det Q{0,0) = 1, or | F p — | 
G p= 1. The formulas (|2l|) are valid for any E, A0 real, k complex, we obtain: 



I oF, (i + e.; x^) ^ -^y^ I oF, + x^) 1 (35) 

to be compared with the similar formulas obtained in ([^) for the absolute values of the 
generalized hypergeometric functions n-^n- 
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In fact, with little changes, one could compute U for any general matrix V with diagonal 
terms equal to zero and non-diagonal terms of general exponential type not neccesarily 
equal. In the particular case of V hermitic: V12 = V21 = k, Uj-ed is unitary and of the form 

Uredizo) = ( ^* -^*] (36) 



with F, G given by Eq.(p^). 

By the unitarity of the matrix (SC), | F p + | G p= 1. And we get the formula: 



17 ^ ^ |2 , I 7 t \ \2 2cosh7re 

I J-^-+J^) I + I J^-+eiix) I = (37) 

For e = this formula reduces to the special case involving the well known Bessel functions 
of order 1/2: J_i/2(a^) = ^2/(7rx) cos(x); Ji/2{x) = yj2/{'Kx) sin(x). 
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Figure 1: The asymmetry A as a function of E/ \m\. Continues line: \ m \ 5 = 2, dashed 
hnes: respectively \ m \ S = 1,1/2, 1/10. Prom A to C: A9 = tt/2, tt, 27r, 67r. For all figures 
^0 = 0. 
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Figure 2: Dependence with the initial angle 9q. Top figure: as in Fig.(|l]) (C) but here for 
^0 = TT- Bottom figure: the quantity A as a function of ^o- = ^tt, E/m = A6/{2m6) 
(peaks in Fig.(|l])(C)). 
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Figure 3: Asymmetry ^ as a function of /S.9 and different \ m \ 5 as before, = 0; 
E5 = 2m5. 
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